Spectral density of an interacting dot coupled indirectly to conducting leads 
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We study the spectral density of electrons p,i CT (aj) in an interacting quantum dot (QD) with a 
hybridization A to a non- interacting QD, which in turn is coupled to a non-interacting conduction 
band. The system corresponds to an impurity Anderson model in which the conduction band has 
a Lorentzian density of states of width A2. We solved the model using perturbation theory in 
the Coulomb repulsion U (PTU) up to second order and a slave-boson mean-field approximation 
(SBMFA). The PTU works surprisingly well near the exactly solvable limit A2 — > 0. For fixed 
U and large enough A or small enough A2, the Kondo peak in pdaiuj) splits into two peaks. This 
splitting can be understood in terms of weakly interacting quasiparticles. Before the splitting takes 
place the universal properties of the model in the Kondo regime are lost. Using the SBMFA, simple 
analytical expressions for the occurrence of split peaks are obtained. For small or moderate A2, the 
side bands of pda(u) have the form of narrow resonances, that were missed in previous studies using 
the numerical renormalization group. This technique also has shortcomings for describing properly 
the split Kondo peaks. As the temperature is increased, the intensity of the split Kondo peaks 
decreases, but it is not completely suppressed at high temperatures. 

PACS numbers: 72.15.Qm, 73.23.-b, 73.63.Kv 



I. INTRODUCTION 

The quantum impurity models, like the Kondo and An- 
derson ones have attracted the solid state physics commu- 
nity due to their complex and rich behavior as well as due 
to their applications to many strongly correlated physical 
systems^ In recent years, the attention into these mod- 
els has increased due to the advances in nanotechnology, 
which for example made it possible to build ideal " single 
impurity" systems with one quantum dot (QD) in which 
the Kondo physics was clearly displayed 2 -'"'^ confirming 
predictions based on the Anderson modeli^i 7 ^. In other 
fascinating experiments, quantum corrals assembled by 
depositing a close line of atoms or molecules on Cu or 
noble metal (111) surfaces have been used to "project" 
the Kondo effect to a remote placed The Anderson model 
used to describe these systems has the additional compli- 
cation of the particular structure of the non-interacting 
states, which cannot be described by a constant density 
of states^ However, it has been shown that the non- 
interacting Green's function can be written as a discrete 
sum of simple fractionsi 11 ' 12 

More recently, systems of several impurities or QD's 
have become a subject of great interest. For example 
non-trivial results for the spectral density were observed 
when three Cr atoms are placed on the (111) surface 
of Au. 13 ' 14 Systems of two^^. three^^ and more 2 ^ 
QD's have been assembled to study the effects of interdot 
hopping on the Kondo effect, and other physical prop- 
erties driven by strong correlations. Particular systems 
have been proposed theoretically as realizations of the 
two-channel Kondo model ) 21 ! 22 the so called ionic Hub- 
bard model^ 3 . and the double exchange mechanism. 24 
In the last years, transport through arrays of a few 
qq^^^ and spin qubits in double QD's 3 ^ have 

been studied theoretically. 



Last year, Dias da Silva et alM proposed a system with 
one interacting QD (1) with a hopping term A to another 
non-interacting QD (2), which in turn is hybridized with 
conducting leads. This hybridization is in principle larger 
than the charging energy of QD 2, in such a way that the 
latter can be considered as non interacting. The model 
is equivalent to an impurity Anderson model in which 
the conduction band has a Lorentzian density of states 
of width A2 . This is one of the simplest variations of the 
model in which the effects of a non trivial structure of 
the conduction band can be studied. It can also be re- 
garded as the simplification of the Anderson model used 
to describe the quantum mirag o 10 ! 11 ! 12 ! 32 or a tight bind- 
ing model for a ring with a quantum do t 10 i 33 when only 
the resonance at the Fermi energy is included. This sim- 
plification is qualitatively valid when the separation in 
energy between resonances is larger than their width. In 
the case of a constant density of conduction states and 
on-site energy Ed below the Fermi energy ep, for large 
enough Coulomb repulsion U, the spectral density of the 
impurity shows two broad peaks (called charge transfer- 
peaks) at Ed < ep and Ed + U > ep, and another peak 
near ep (the Kondo peak). 

Previous calculations of the spectral density of the im- 
purity in the model for the quantum mirag o 11 ' 32 ^ and 
in confined structures 3 ^ have shown that the Kondo peak 
near the Fermi energy splits in two for certain parame- 
ters. Dias da Silva et al. studied the spectral density 
at the QD 1 pd a (uj) and magnetic susceptibility of the 
symmetric model using numerical renormalization group 
(NRG). In agreement with previous studies, they find 
that when A increases beyond a certain critical value A c , 
the Kondo resonance splits. In a similar way, for large 
enough A split peaks are observed if U is smaller than 
a critical value U c . In fact previous calculations for the 
spectral density of an impurity inside a circular corral 
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show split peaks in the non-interacting case U = that 
turn into one peak as U is increased^ The splitting con- 
dition is proposed to correspond approximately to the 
equality \/2Tk — A2, where Tk is the Kondo temper- 
ature. However, the results presented to support this 
reasonable statement are limited. The spectral density is 
shown only for three particular sets of parameters. More- 
over, the ability of the NRG to describe peaks far from 
the Fermi energy might be questionable because due to 
the logarithmic discretization, the resolution at a given 
energy u> is proportional to ui — ep. In addition, the spec- 
tral broadening affects the resolution^ 

In this paper we use perturbation theory in the 
Coulomb repulsion U (PTU) up to second order— ^ and 
a slave-boson mean-field approximation (SBMFA)^ to 
study in more detail the spectral density, the conditions 
for the appearance of split peaks near ep, and the evo- 
lution of the energy scale Tk with the parameters of the 
model. We also compare with exact results in the limit 
A2 — > 0, and derived from the Friedel sum rule4£ As 
explained in the next section, both approximations take 
a simpler form in the symmetric Anderson model con- 
sidered, and the PTU works better in this case, allowing 
us to describe new physics, and obtain reasonably robust 
results with comparatively simple mathematics. 

For the SBMFA we used a straightforward extension 
of the method proposed by Kotliar and Ruckenstein to 
the Hubbard models 

The paper is organized as follows. The model, the ap- 
proximations (PTU and SBMFA) and some exact results 
derived from Fermi liquid properties and the case A2 = 
are presented in Section II. Several figures illustrating the 
main results are presented in Section III. Section IV is a 
short summary and discussion. 

II. HAMILTONIAN, APPROXIMATIONS AND 
EXACT RESULTS 



B v = J2Wi(aUi. + tt-c). (2) 

la 

The sum H 2 + Hi + Hy is a non-interacting Hamiltonian 
which can be put in the diagonal form J2ka e k c ka Ckcr by 
means of a canonical transformation. In this basis the 
Hamiltonian takes the form of the impurity Anderson 
model for a general band structure and hybridization 

H = H 1 +J2 ^c{ a c k a + Vk(4c k<7 + H.c), (3) 

k<7 k<7 

where V k = Aja^, c\ a }. 

The spectral density of electrons at the interacting QD 

is 



Pda M = — [Gdda (W - if]) - Gdda (w + IT))] , (4) 

where 77 is a positive infinitesimal and calling z = lo + irj 
(z = u) — ir)), the retarded (advanced) Green's function 
at the interacting QD can be written in the for m 11 ' 40 



<^dda\Z — ^ ^ i—s r; TTi loj 

Z — Mji — l^ s da{z) - hdda(Z) 

where 'Sdda(z) is the self-energy due to the interaction U 
and 



Z sd a(z)=Y<—- W 

k 

While this result holds for a general Hamiltonian of the 
form J3J, for the particular case of Eqs. ([I]) and using 
equations of motion,— it can be shown thai^ 



A. Model and equations for the spectral density 

The system is described by the following Hamiltonian: 



H = H 1 +H 2 + H l +H x + H v , 



(1) 



where Hi {H2) describes the interacting (non- 
interacting) QD, Hi the leads, and the last two 
terms are the hybridization of the non-interacting QD 
with the interacting one and the leads respectively. In 
standard notation 

H l = E d ^dld a + Ud\di i d\d i , 

(7 

Hi = e a a\ag-, Hi = eicj a cia, 

a la 

H x = \J2(dUa +H.c), 



SWz) = A 2 G° Qaa (z), 



(7) 



where G aaa (z) is the Green's function of the non- 
interacting QD in the absence of the interacting one (for 
a Hamiltonian H 2 + Hi + Hy). The model assumes con- 
stant values for the matrix element Wi and the density 
of states of the leads These assumptions are usu- 
ally very good approximations for the range of energies 
of interest in QD's. Calling A2 = 7r/9|IU| 2 , they lead to 



uj - e a ± iA 2 



(8) 



In the following, as in Ref. [31], we take the origin of 
one-particle energies at the Fermi energy {ep = 0), and 
take e a — 0, = —U/2, corresponding to the symmetric 
Anderson model. In summary, the model is equivalent 
to an impurity Anderson model with the hybridization 
function 
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A(uj) = -ImS s d (7 (w + irf) = 



uj 2 + A 2 ' 



(9) 



Note that in general, the real part of T, sda can be ob- 
tained from the relation 



A(w) = ReE sdrT (iu) = P 



A{e)de 
tt(u> — e) 



(10) 



If some approximation for Tt dda {z) is used, the above 
equations (OJ , (J5J) , ([TJ) and (J8J) define the spectral density. 



B. Approximations 

1. Perturbations in U 

The starting point of the perturbation theory in U 
(PTU) is a non-interacting problem which includes some 
one-body potential so that the effective on-site energy 
of the interacting electrons E d is modified to E e JJ . A 
possible choice is the Hartree-Fock value E e JJ — Ed + 
U{dUa}- This one-body potential is compensated in the 
perturbation, which includes it (with the opposite sign) 
in addition to the interaction term £/djdfdjdj. In gen- 
eral, there are better choices for E e J^ which have been 
used recently in several calculations of transport prop- 
erties in nanoscopic systems, including non-equilibrium 
situations and applied magnetic fields.— However, for 
the symmetric Anderson model at equilibrium without 
applied magnetic field (for which {d\d a ) = 1/2) these 
choices coincide with the Hartree-Fock result and E e JJ = 
Ed + U/2 = ep. In this case, for a flat band the theory 
is quantitatively correct up to U ~ 8A, where A is the 
resonant level width i 42 ' 43 Working up to order U 2 we can 



[G dd a(z)}- 1 = [ 9dda (z)}- 1 -J: d 2 l(z) ) 



(11) 



where g dda (z) is the Green's function for the interaction 
treated in Hartree-Fock (including contributions of first 
order in U to S^dcr) 



9dda{z) 



Z - Y, sd a{z) 



(12) 



(2) 

and Y/ dda {z) is the contribution to the self-energy of sec- 
ond order in U evaluated from a Feynman diagram in- 
volving the analytical extension of the time ordered gdda 
to Matsubara frequencies 

^dLi iuJ h T ) = U 2 T^2 9dda(iui - w m )x;(w m ),(13) 

m 

Xo(iv m ) = -T22g dda (iu> n )g Ma (iu) n + iis m Xl4) 



where u> n = (2n + l)irT and v m = 2rmrT. 

Our task is simplified because \<r can t> e calculated 
analytically due to the fact that g dda can be expressed as 
a sum of two simple fractions 



gdd<r{u+iri) = 



uj + iA 2 



lu 2 + iA-2LU — X 2 ^—^ lu — ri + iSi 



ai 



(15) 



This is the retarded Green's function. In order to evalu- 
ate perturbative diagrams, like Eqs. (fT3| . (fT4|l . one needs 
the analytical extension of the time ordered Green func- 
tion to imaginary frequencies u) — ► io ^ 11 ' 44 : 



gdd, 



ai + ai + (a; - a;)sgn(u„) 
2 [w n + i^sgn(w n ) - n] 



(16) 



where ai is the complex conjugate of ai and sgn(x) is the 
sign of x. 

The sums over Matsubara frequencies, Eqs. (fT^)) . (fTi)) 
have been done in Rcf. The latter can be expressed 
in terms of the digamma function ^(^ii The final result 
for the retarded quantities is 

E ddl( w ) = — d y [ y/r „ i 9dd*(u - y)lmxcr(y) + 



1 



ev l T + 1 



Img dd(T (-y) x<t(u + y)], (17) 



where 

/•7T ' ■ 



I in 



- (Bi m (uj) + A hn (-uj))^i(uj) 

(18) 



with 



^im(^) = 
Bl m {i0) = 



ai a r , 



w + n - r m - i<5i + £<5„ 

to - n + r m + iSi + iSr, 
1 Si + i(n - u) 



2-kT 



(19) 



Eqs. Q, Q, ©, Ijnjl, ((UJ) and fTfjl to (fl9l define the 
spectral density of the interacting dot within PTU. 



2. Slave bosons 

The basic idea of the slave boson formalism of Kotliar 
and Ruckenstein 3 ^ is to enlarge the Fock space to in- 
clude bosonic states which correspond to each state in 
the fermionic description at the interacting QD. The vac- 
uum state at this site is represented as e> |0), where 
et is a bosonic creation operator corresponding to the 
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empty QD; similarly syf£ |0) represents the singly oc- 
cupied state S a |0), where s\ is a bosonic operator for 
singly occupied sites with spin a, and fl is a fermion op- 
erator. The doubly occupied site djdj |0) is represented 
as b^ftf] |0) In this way the interaction term can be ex- 
pressed in terms of boson operators as Ub'b and the in- 
teractions between fermions disappear from the Hamil- 
tonian. In other words, the fermion operator d\ j s ex _ 
pressed in terms of fermion operators that do not inter- 
act between them as d\ = (tfs s + s£e)/j!. The bosonic 
operators should satisfy the following constraints 



e f e + 



^4s CT + 6+6 = 1 



(20) 



Introducing Lagrange multipliers for these constraints, 
the Hamiltonian takes the form 



[R*{b ] s s + 4e)L ff /ta CT + He] 



-n{Je + ^2 4 s <r + tfb - 1) 



(21) 



Here, the factors 



(1 



e t e _ s t SCT) -i/2 



tfb) 



-1/2 



(22) 



are equivalent to one when treated exactly, but they are 
introduced in such a way that in the slave-boson mean- 
field approximation (SBMFA) the correct result in the 
non-interacting case is reproduced. 

In the SBMFA, all the boson operators are replaced by 
numbers and the free energy is minimized with respect 
to them. In the symmetric Anderson model, without ap- 
plied magnetic field the problem simplifies considerably. 
In this case s 2 is independent of spin and e 2 = b 2 . Also 
(d\d a ) = {flfo) = 1/2. From here E d -fi a = e F = and 
using Eqs. (j20f s 2 = 1/2 — b 2 and only one independent 
variable remains. Also R a = L a = \/2. The change in 
free energy due to the impurity can be written a*^ 



-Im 



f(u)lnG ffa (uj + irj) du; + Ub 2 

+2 A v(s 2 + 6 2 ),(23) 



where f(u>) is the Fermi function and Gff a (u>) is the 
Green's function of the / operators for any spin. The 



Hamiltonian takes the same form as that of a non- 
interacting problem with Ed = ep — and renormal- 
ized hopping A = R a (b^ s s + s\e)L a \, that in our case 
simplifies to 



A = 46\/l/2 - 6 2 A. 
Then, as in the previous subsection 



G ffa (u + iri) 



iA 2 



iA 2 uj - X 2 



(24) 



(25) 



Decomposing this expression in simple fractions [as in 
Eq. (fTSf ] and replacing in Eq. (j23|) , the integral can be 
evaluated analytically at zero temperature. The result 
should be separated in two cases depending on the sign 
of 2A — A2. Except for irrelevant constants, the result is 



AF 



-{ln(l + r 2 ) + r(2 arctanr -1 - tt)} + Ub 2 



for x > 1, and 

AF = — {(l-r)ln(l-r) + (l + r)ln(l + r)} + [/6 2 
for x < 1, with 



r = yj\x 2 - 1|, x = 2A/A 2 



(26) 



Minimizing AF defined by Eqs. (J26J) and (|24[). one 
obtains a transcendental equation for b. After solving 
this, a characteristic energy scale or Kondo temperature 
can be defined by the gain in energy with respect to the 
unhybridized case: 



(27) 



Tk = AF(A m in) — AF(0), 



where A m i n is the value of A evaluated with Eq. (|24|) for 
the value of b that minimizes the energy. 

The spectral density for real d electrons near the Fermi 
energy becomes, using d\ = i? cr (6^s 5 - + sJ.e)i CT /J in the 
SBMFA 



G dda {u + w?) ~ (A min /A) 2 G //CT (w + in). (28) 

In the resulting spectral density pdaiyj) [calculated from 
Eq. dl])], the "charge transfer" peaks near E d and E d + U 
are lost, as explained in more detail in the next section. 

From the change of sign of d 2 p da (oj) / doj 2 evaluated 
at the Fermi energy, one finds that (in a similar way as 
in the non-interacting case^) the critical condition for 
the appearance or disappearance of split peaks near the 
Fermi energy is 



V2A„ 



(29) 



Combining this equation with Eq. (|24p and the min- 
imization condition dAF/db — 0, we find in order to 
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have split peaks, it is necessary that A > A2/V2 and 
that U <U C , where 



8A^/ A| 
A 2 I 2A 2 



1/2 



(30) 



This analytical result can be inverted to give the mini- 
mum value of A required to have split peaks for fixed U 
and A 2 : 



A, = 



U 



u 2 

4A^ 



1/2 



1/2- 



(31) 



C. Exact results 



1. Fermi liquid properties 



d\ — > d a , cl a — > -c k < CT , with ek' = -£k- Calling 
the thermodynamic potential and |n) a complete set of 
eigenstates, we can write^ 



m ) 



; e -0E n+e -fSE m 
LO + E n - E m + %Tj 

(36) 



Because of symmetry one has (H^<t 



(Tn\TdtT^\Tm) = (Tn\d a \Tm) = (Tm\dl\Tn), 
and the eigenstates \m) and \Tm) have the same energy. 
Then changing the labels of the sum above \n) by \Tm) 
and |m) by |Tn) 



G dda {-u-ir,) = e^ n Y,\(Tm\dl\Tn)\ 



e -0E n + e -p Em 
-LJ + E m - E n - it] ' 



and using the symmetry property of the matrix elements 
indicated above, one finds 



It is known that at the Fermi energy, the imaginary 
part of the self energy due to interaction vanishes in a 
Fermi liquid^ Using ^^^(ej?) = 0, and Eqs. ^ and 
([5]) with Y^ sda (uj ± irj) — A a (oj) =F «A ct (oj), the spectral 
density at the Fermi energy can be written in the form 



COS lf a 

7rA CT (e F )' 



(32) 



where 



arctan 



( E d + T, ddr7 (e F ) + A tT (e F ) 



A a (e F ) 



(33) 



In addition, for the general Anderson model, the Friedel 
sum rule is valid^ 



= 7T ((did,) - ~ 



9o 



-lm / duiG dda (uj + IV)- 
J — 00 



du> 



(34) 



As explained below, for the symmetric Anderson model 
without applied magnetic field, with the choice ep = 0, it 
can be shown that the real (imaginary) part of G dda {ui + 
irj) is odd (even). Then Eq. ^ implies Re [£y+£ s d CT (0)+ 
E dd(7 (0)] = 0. From Eq. <p a = and from Eq. 



1 



ttA(O)' 



(35) 



The symmetry properties of G dda {uj + iff) can be 
demonstrated using the Lehman representation of the 
Green's function^ and the fact that the symmetric An- 
derson model is invariant under the transformation T : 



G dda (-uj - irj) = -G d dcr{u + irj). 



(37) 



We have verified that Eqs. (fM]) and (|33|) are satis- 
fied by the approximations (PTU and SBMFA) presented 
above. For the PTU the second member of Eq. (|3"4"1) has 
been evaluated numerically and we find that it is zero 
within the accuracy of the computer. Instead, within 
NRG, Eq. ([55]) is satisfied only approximately^ 



2. The limit A 2 







This limit coincides with the so called atomic limit of 
the Anderson model, reported previously in Appendix C 
of Ref. Q] and in Ref. WA We describe here the main 
results for the symmetric Anderson model. 

The ground state is a two-particle singlet with energy 
Eg = E t — T K , where E t = —U/2 is the energy of the 
triplet state and the characteristic energy Tk is 



T - U 




(38) 



This energy coincides with the gain in energy due to 
hybridization, Eq. (f27f . Note that for U — > 00, 
Tk = 8X 2 /U in contrast to the result for a flat band 
T K ~ Dexp[-7rZ7/4A]. 

The spectral density is easily calculated using Eq. (g]) 
and the Lehman representation Eq. (j36|) of the Green's 
function. For later use, we display here the result at 
T = 



p da (u) = A[5(uj - E A ) + 5{u + E A )] 

+ (1/2 - A)[8(w - E B ) + 5(lu + E B )],(39) 



6 



where ±Ea are the positions of the peaks nearest to the 
Fermi level and A their weight, while ±Eg and 1/2 — 
A are the position and weights of the "charge transfer" 
peaks nearer to ±£7/2. The values of these quantities are 



E A = 



U 



U 



E B = - 
A = - 




y/(l + y 4 )(l - m) + V(i - wKi + 2/s) 

21-1/2 



W = [l + (jA/C/) 
For f7 — * oo, they simplify to 



(40) 



to the extra broadening of these peaks was presented by 
Logan et al42. These authors also show that qualitatively 
these peaks can be understood using an alloy analog ap- 
proach (AAA) that consists in replacing the system by 
a non interacting system in which the on-site energy of 
the localized electrons is either Ed or Ed + U with equal 
probability. This approach misses the Kondo peak but 
describes surprisingly well the charge transfer peaks for 
large In our case, the AAA also shows high and 

narrow peaks for small enough A2 and the reason is that 
A(cj) is considerably smaller for uj — Ed or u> = Ed + U 
than at the Fermi energy A(0). While this analysis gives 
confidence to our results, the reader might still wonder 
if the aspect of the side peaks is due to a shortcoming 
of the PTU. For constant A(w) (large A2) the approach 
is quantitatively valid for U < 8Aj 42 i 43 In the following 
subsection we analyze the precision of the method for 
small A2. 



E A = 6X 2 /U 
E B = U/2 + 10X 2 /U 
A = 18\ 2 /U 2 . (41) 

When the temperature is increased, more peaks appear 
in the spectral density and the weight A of the peaks 
nearer to the Fermi energy decreases. 

III. APPROXIMATE RESULTS 

A. General aspect of the spectral density 

In Fig. [T]we show the spectral density calculated with 
PTU. We have fixed A(0) = A 2 /A 2 as the unit of energy. 
Then according to Eq. ([35)1 the spectral density at the 
Fermi level is 1/tt, for all sets of parameters. On can 
see that the PTU satisfies this condition derived from 
Fermi liquid properties. In addition to this observation, 
there are two other noticeable properties that can be ob- 
served in the figure, in comparison with the case of con- 
stant A(w). One is the presence of split peaks near the 
Fermi level for large enough A and U, already observed 
befor o 12 ' 31 as discussed in the introduction. The other is 
that the "charge transfer" peaks at Ed and Ed + U are 
unusually high and narrow, particularly when split peaks 
appear. We begin discussing the latter fact. This has not 
been noticed by previous NRG calculations. 31 We believe 
that this is due to the lack of resolution of the NRG for 
structures that are far form the Fermi energy. We re- 
mind the reader that for constant A(cu) these side peaks 
are broad and of smaller amplitude than the only cen- 
tral peak (the so called Kondo peak). In particular, in 
the symmetric case Ed — —U/2, while the width of the 
Kondo peak (the half width at half maximum) is of the 
order of the Kondo scale Tjc and its height is 1/(ttA), the 
width of the side peaks is of order 2A and their height 
is near l/(47rA). A discussion of the effects that lead 




oVA(O) 



FIG. 1: Spectral density of states fo r U — 4 and several values 
of A 2 /A(0) keeping A = x /A 2 A(0). 

Near the Fermi energy, we see that as the ratio A/A2 
increases, first the Kondo peak narrows, then its splits 
in two and with further increase in A/A2 the split peaks 
became narrower and higher. They also tend to move 
towards the Fermi energy but they never merge into one 
again (for A » A 2 , the exact results of Section II. 2 be- 
come qualitatively valid). 

As expected, if U is lowered, the values of A/A 2 needed 
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FIG. 2: Spectral density of states for U 
of A 2 /A(0) keeping A = v /A 2 A(0). 



2 and several values 



to have split peaks are smaller. This can be seen compar- 
ing Fig. [T] (U = 4 A) with Fig. CU (f7 = 2A), or Fig. 2 of a 
previous work on the quantum mirage effect^ in which 
split peaks are obtained only in the non-interacting case 
of the model used. 
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B. Variation with A 2 

In Figs. [3] and S] we show the evolution of the spectral 
density calculated with PTU as A2 is increased, start- 
ing from very small values. As A2 increases, all peaks 
broaden and those nearer to the Fermi energy approach 
each other, merging into one for large enough A2. This 
evolution is similar to that reported previously in Fig. 2 
of Ref. as a function of the width of the resonances 
in the Anderson model for an impurity inside a quantum 
corral. In Fig. [3]we choose A as the unit of energy and the 
evolution ends with only one peak at the Fermi energy 
when A2 reaches 1. Increasing A to 1.5, the magnitude 
of the splitting increases and the peaks near the Fermi 
energy are higher and remain split for A2 = 1. 

For A2 = 0.02 we have integrated numerically the spec- 
tral density pda{<+>) below each peak to compare with the 
exact results given by Eqs. (|40[) . From the energy of the 
maxima of pd<j{w) (±Ea, ±-Eb) and the resulting values 
of the peak weights (A and 1/2—^4) we obtain Ea = 0.82, 
Eb = 3.65, and A — 0.325 for the parameters of Fig. [3] 
and E A = 1.36, E B = 4.96, and A = 0.401 for the pa- 
rameters of Fig. [4j The corresponding exact values for 
A 2 = given by Eqs. (gO]) differ in less than 0.004 from 
those given above. This shows that the PTU works very 
well for small A2 and provides confidence to the results 
presented above. 

In Fig. 0(a), the result of the SBMFA is also shown for 
comparison. This approximation also gives qualitatively 
correct results for the two peaks nearer to the Fermi en- 



FIG. 3: Spectral density of states for U = 4, A = 1 and several 
values of A2. 

ergy for small A2, but quantitatively the PTU is superior 
in this limit. 

In Fig. we show the effect of a larger U. When 
comparing with Fig. [H we see that the peaks near the 
Fermi energy become closer again and spectral weight is 
transferred to the side peaks. In this case, we obtain 
within PTU for A 2 = 0.02, E A = 1.09, E B = 6.10, and 
A = 0.261 . These values differ from the exact results for 
A2 = in less than 1%. 



C. Description of the splitting in terms of 
quasiparticles. 

While the limit A2 — > provides a qualitative descrip- 
tion of split peaks near the Fermi energy, a more precise 
description can be given in terms of weakly interacting 
quasiparticles, for which the analytical results of a non- 
interacting systems with renormalized parameters gives 
an accurate description as a starting point, as we show. 

Hewson has noticed that near the Fermi energy, it is 
possible to reformulate the PTU in terms of quasipar- 
ticles, for which the residual interactions are smaller^ 
The starting point of renormalized perturbation theory, 
as in Fermi liquid theory, is an expansion of the self en- 
ergy Sddo-(w + ii]) around the Fermi energy (zero in our 
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case). Since ImSdd cr (w + irj) 



near the Fermi energy 



[see Fig. [6] (b)] , to linear order in u> we can approximate 
Eq. © as 



1 



w(l - dReEdto/duj) -E' d - S sdCT (w) 



w - zE' d - zE sd(T (a;) ' 
where E' d = Ed + Srfdo-(O) (zero in our case) and 

1 



(42) 



(1 - dRei: dd(7 /dLu) 



(43) 



The interactions however still modify this picture and a 
renormalized self energy should be introduced. 49 How- 
ever, we find that Eq. (142(1 with z calculated numerically 
from the slope of ReS^do- leads already a very good de- 
scription of the spectral density near the Fermi energy. 
This is shown in Fig. [S] for two different sets of param- 
eters, one leading to split peaks and the other not. The 
discrepancies between the complete PTU results and this 
approximation are appreciable only sufficiently far from 
the Fermi energy, where the influence of the side peaks, 



ignored in Eqs. (|42|) cannot be neglected. The result- 
ing values of z for A2/A(0) from 0.2 to 1 in steps of 0.2 
with A = v /A 2 A(0) are 0.20, 0.28, 0.33, 0.36 and 0.38 
respectively. 

Note that this approximation is equivalent to use a 
non-interacting model with renormalized hybridization 
A = s/zX, [sec Eqs. (J9j) and (fit))) ] and decrease the re- 
sulting Gddoiu + if]) by a factor z < 1. This is very 
similar to the SBMFA [see Eqs. fl23]) and (|2"g|)]. except 
for the fact that in the latter, the renormalized value of 
A is obtained by a minimization of the energy and not 
by deriving an approximation to the self energy. Thus, 
the results presented above can be regarded as a support 
to the qualitative validity of the SBMFA to describe the 
spectral density pdo-(tjj) near the Fermi energy. A com- 
parison between results of PTU and SBMFA for pda(<jj) 
is presented in subsection F. 



D. Critical conditions for split peaks 

Within the SBMFA we have obtained analytical results 
[see Eqs. (|30l) and (|3"Tj) ] for relations between the param- 
eters when split peaks just appear. This leads to a kind 
of "phase diagram" for the appearance of split peaks. 
The boundary line is shown in Fig. [7J Split peaks are 
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FIG. 6: (a) Spectral density of states near the Fermi level for 
two values of A2 keeping A(0) = 1. Full lines: [7 = 4, dashed 
lines: renormalized non-interacting case, (b) Self energy for 
A 2 = 0.2 
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FIG. 7: Boundary separating the region of parameters for 
which split peaks in the spectral density are present. 
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expected for larger values of A or smaller values of U , the 
upper left of the curve. This results should be regarded 
as qualitative. Comparison with Figs. [T] and [2] show that 
slightly larger A c are expected within PTU. In fact, while 
using the factors R a and L a [see Eqs. ([22])] in the SBMFA 
leads to the correct results in the non-interacting limit 
U ^ 0, these roots should be eliminated to obtain the 
right exponential dependence of Tk in the limit U — > 00 
for constant A(w). 2 - This modification increases A c by a 
factor 2. Therefore, we expect that the actual value of 
A c is larger than that indicated in Fig. [3 particularly for 
larger values of U. 

It has been suggested that in general, the critical con- 
dition can be defined as \[2T K ~ A2, where T K is the 
Kondo temperature at the critical line^i The first mem- 
ber of this equation as a function of U is represented in 
Fig. [HJ with the Kondo temperature defined by Eq. ([2"?)) . 
Within a factor 3 we find that this condition is correct, 
which is not bad in view of the exponential dependence 
expected for Tk on the parameters near the transition. 
We note that due to a cancellation of factors, the same 
T K is obtained if the roots R a and L a are dropped in the 
SBMFA. Therefore, this result seems to be robust. 



FIG. 8: Kondo temperature as a function of U at the bound- 
ary for the appearance of split peaks. 



E. Kondo temperature and universal behavior 

For constant A(w) in the symmetric Anderson 
model, the Kondo temperature is given by Tk — 
D exp[— 7rE//4A], where D is the band width, and the 
properties of the system, depend on Tk rather than U 
and A individually. Thus, the behavior is universal in the 
sense that different curves can be mapped into one with 
appropriate scaling. In the present case, when A 2 3> Tk, 
A(oj) is flat on the scale of Tk, and one expects the same 
behavior with A replaced by A(0) = A 2 /A2. Therefore, 
in this region, when in general there are no split peaks, 
one expects Tk — -Dexp[— 7rJ7A2/(2A) 2 ]. However, for 
A2 — * 0, one has a different behavior, given by Eq. (|38|) . 

In Fig. [5] we show ln(TRr/Aa) as a function of the 
ratio UA2/X 2 for several values of A, calculated with the 
SBMFA. As expected, for large values of this ratio, for 
which there is only one peak near the Fermi energy, all 
curves merge into one straight line, indicating a universal 
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FIG. 9: Natural logarithm of the Kondo temperature as a 
function of the ratio UA2/X 2 for several values of A/ A2 . From 
bottom to top 0.25, 0.35, 0.5, 1, 1.5, 2, and 2.5 



behavior, and an exponential dependence of Tk on this 
ratio. However, as UA2/X 2 decreases, already for rather 
large values (~ 10) for which no split peaks are expected 
according to Fig. [7] the values of Tk for different A 
start to deviate between them and from the exponential 
behavior. The decay is faster (slower) than exponential 
for the larger (smaller) values of A considered. 



F. Comparison of the spectral density obtained by 
different methods 

In Fig. [10] we show the spectral density calculated by 
PTU for the same parameters as those used in recent 
NRG calculations^ The scale is chosen to display the 
high intensity of the side peaks. 

In general, it is difficult to determine the range of valid- 
ity of the PTU in terms of the parameters U, A and A2. 
For A 2 3> Tk (well inside the region without split peaks) , 
comparison with results of quantum Monte Carlo*2 and 
other calculations^ suggest that PTU is quantitatively 
valid for U < 8A0, where Ao = A 2 / 'A2. This condition 
suggests that PTU is near its limit of validity for the 
smallest value A = 0.0354 considered in Ref. [3l|. When 
A 2 < Tk, one might expect that the above condition 
should be replaced by U < 8A av , where A a „ < A is 
A(w) [see Eq. © ] averaged over a range of energies of 
the order of Tk- In any case, comparison with the results 
of the SBMFA discussed below for the other two values 
of A used, indicates that the results of the PTU are reli- 
able for these sets of parameters. The comparison with 
exact results for A 2 <C A discussed in subsection B also 
supports the validity of the PTU at zero temperature. 

For A = 0.0354, Pda(u) reaches 90 in the arbitrary 
units chosen in Ref. l3ll . These narrow peaks near —{7/2 
and U/2 are probably lost by the NRG approach, which 
shows only broad features there. However, these peaks 



might be experimentally accessible and are therefore rel- 
evant. 
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FIG. 10: Spectral density of states for U = 0.5, A2 
and several values of A calculated with PTU. 



0.02 



The density for the larger value of A is shown in Fig. 1111 
together with the corresponding SBMFA result. The side 
peaks are lost in the latter approximation. Both results 
agree qualitatively. We believe that as discussed in sub- 
section D, the use of the roots R a and L CT [see Eqs. ([22])] 
in the SBMFA leads to renormalization factor z = (A/A) 2 
which is larger than the correct one for large values of U. 
This is probably the reason of the discrepancy. In any 
case both approaches lead to peaks around the Fermi 
energy which are three or four times larger than those 
reported in Ref. [3l|. This again is likely due to lack of 
resolution of the NRG results for energies away from the 
Fermi energy, as argued below. Note that reducing A 2 
to zero, the exact solution [see Eqs. ([22]) ] gives delta 
functions at u> = ±0.053 and u> — ±0.347, with weights 
A = 0.221, and 1/2 - A = 0.279 respectively. This re- 
sult is quite consistent with the PTU one shown in Fig. 
[TT1 For example, a lorentzian fit of a side peak obtained 
with PTU gives a center at ±0.343, an area 0.284, and 
a total width 0.0147 (smaller but of the order of A 2 ). 
Furthermore, taking these four delta functions artificially 
broadened as it is usual in NRG calculations (using Eq. 
(10) of Ref. Hil with b — 0.6) and convoluting it with 
a Lorentzian of total width A 2 we obtain the result dis- 



which is very similar 
3ll with central split 



played by a dotted line in in Fig. QT 
to that shown in Fig. 2 (c) of Ref. 
peaks of intensity near 2 and only very broad bumps re- 
placing the side peaks near oj = ±0.35. These arguments 
indicate that the NRG, at least in its usual form, is in- 
adequate to describe the spectral density in the region of 
parameters in which split peaks are present. 

The results for the spectral density within PTU and 
SBMFA near the Fermi energy for the smaller two values 
of A are shown in Fig. [T2] For the same reason discussed 
above, we believe that the splitting and intensity of the 
peaks for A = 0.0627 within the SBMFA are exaggerated. 
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FIG. 11: Spectral density of states near the Fermi energy for 
U = 0.5, A 2 = 0.02 and A = 0.078. Full line: PTU, dashed 
line: SBMFA, dotted line: result for A2 = broadened as in 
NRG calculations (see text) 



As before, the PTU results predict narrower peaks near 
the Fermi energy than the NRG results. For the smallest 
value A = 0.0354, we believe that both, the PTU and 
the SBMFA including the roots R a and L a exaggerate 
the width of the only peak near the Fermi energy and 
the NRG is expected to provide more reliable results. In 
comparison with the latter, the SBMFA without the roots 
seems to predict a too narrow Kondo peak. However, 
even in this case, the NRG is not able to describe properly 
the narrow side peaks. We expect that this defficency 
persists as long as A 2 < U/2. 
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FIG. 12: Spectral density of states near the Fermi energy for 
U = 0.5, A 2 = 0.02 and two values of A. Full line: PTU, 
dashed line: SBMFA. 



G. Effects of temperature 

In Fig. [13] we show the evolution of the spectral den- 
sity with temperature in a region in which there are split 
peaks, calculated with the PTU. Similarly to the case of 
constant A(w), the peaks near the Fermi energy decrease 
in a temperature scale of the order of the half width at 
half maximum of the structure. In contrast to that case 
however, some structure remains even at the highest tem- 
peratures. This can be understood qualitatively from the 
behavior observed in the exactly solvable limit A2 = 0. 
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FIG. 13: Spectral den sity of s tates near the Fermi level for 
A 2 = 0.2A(0), A = y/AzA(Q), U = 4 and several values of 
the temperature T/A(0). 

In addition, the evolution of the side peaks with tem- 
perature is more marked than in the case of constant 
A(w), probably due to the fact that the broadening ef- 
fects of the temperature are more noticeable when the 
peaks are narrower. 

At temperatures higher than T = A(0), there is little 
variation in the spectral density. The result for T — 
5A(0) (not shown) coincides with that for T = A(0) 
within the width of the line in the figure. It is curious 
that some structure appears near w/A(0) = ±1.2 at high 
temperatures. Its origin is unclear to us. 



IV. SUMMARY AND DISCUSSION 

We have studied an impurity Anderson model in the 
symmetric case, in which the hybridization function A(u>) 
has a Lorentzian shape of width A 2 . The model can be 
regarded as a simplified version of that proposed to de- 
scribe the projection of the Kondo effect in quantum cor- 
rals (the quantum mirage)?^ and is relevant for systems 
of double QD's in which one can be regarded as non- 
interacting and connected to metallic leads (or other sys- 
tem of non-interacting electrons) and the other is inter- 
acting and connected as a side dot to the non-interacting 

r,„ a 3i 
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The characteristic low-energy scale or Kondo temper- 
ature Tk of the model, changes from the known ex- 
ponential dependence Tr — D exp[— 7rf7/4A(0)] in the 
Coulomb repulsion U for large U and a flat A(u>) (large 
A 2 ), to T K = 8X 2 /U for large U and A 2 -> 0, where A is 
the interdot hybridization. For small enough A 2 < Tr-, 
the scaling properties of the Kondo regime are lost, but 
by extension we continue to call "Kondo peak" the struc- 
ture near the Fermi energy in the spectral density. 

When A 2 becomes smaller than ~ TV/2 (see Fig. [5]) 
the characteristic Kondo peak in the spectral density 
Pdo-(w) of the interacting QD near the Fermi energy splits 
in two. This splitting was reported before in the context 
of the mirage effect (but not discussed in detail there) 
and in recent calculations using the numerical renormal- 
ization group (NRG) 4i Our calculations using perturba- 
tion theory in U and a slave-boson mean field approxi- 
mation (SBMFA) agree qualitatively with these results, 
but predict narrower split peaks. Moreover, we find that 
the side peaks of pda(u) near the charge transfer energies 
are much narrower and higher that in the usual case of a 
flat A(u>) and than the results presented in Ref. [3l|. 

In Section III F, we have provided several arguments 
that indicate that the above mentioned discrepancies are 
due to shortcomings of the numerical technique. These 
can be partially overcome if density matrix renormal- 
ization group is used in combination with NRG. 26 An- 
other way of improving the NRG results which has been 
shown to lead to sharper charge transfer peaks is to cal- 
culate the self energy as a ratio of two Green's functions 
^dda(z) = U{{d a dldv,dl)) z /{{d a ,dD) z , and replace the 
result in Eq. 36 This has the advantage that the effect 
of Ttsdcriz) is taken into account exactly. 

The splitting of the peaks can be understood quali- 
tatively from the limit A 2 — * 0, and quantitatively in 
terms of weakly interacting renormalized quasiparticlcs, 
as described in Section III D. 



The SBMFA provides an analytical expression, lead- 
ing to a diagram for the region of parameters for which 
split peaks in pd a {u>) are expected. This is represented 
in Fig. [7J In the region in which split peaks are present, 
some structure remains near the Fermi energy even at 
temperatures much higher than the Kondo temperature 
T K - 

The approximations that we have used satisfy Fermi 
liquid relations and work well in the limit of small A 2 . 
In particular the agreement of perturbative calculations 
with the exact results for A 2 — > at zero temperature is 
surprising. 

If the non- interacting resonance (the energy of the non- 
interacting QD) is shifted away from the Fermi energy in 
some energy larger than its width A 2 , we expect a dra- 
matic change in the spectral density at low energies, with 
a very narrow Kondo resonance at the Fermi energy, due 
to the decrease of the non-interacting density of states 
at the Fermi level, and the exponential dependence of 
the Kondo energy scale with this density. This is based 
on previous calculations using approximations^ and sup- 
ported by NRG results^ Instead, the side peaks should 
not be affected substantially. 

The presence of narrow side bands in the region in 
which the Kondo peak splits in two can in principle be 
tested experimentally in QD systems. The spectral den- 
sity at the interacting dot can be measured in transport 
experiments in which another lead is added i 50 ' 51 
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